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Abstract. Acellular slime mould Physarum polycephalum is a monstrously large sin- 
gle cell visible by an unaided eye. It shows sophisticated behavioural traits in foraging 
for nutrients and developing an optimal transport network of protoplasmic tubes span- 
ning sources of nutrients. The slime mould sufficiently approximates man-made trans- 
port networks on a flat substrate. Does slime mould imitate man-made transport net- 
^ works on three-dimensional terrain as well as it does on a flat substrate? We simplified 

Q the problem to approximation of a single transport route. In laboratory experiments 

with 3D Nylon terrains of USA and Germany we imitated development of route 20, 
the longest road in USA, and autobahn 7, the longest national motorway in Europe. 
We found that slime mould builds longer transport routes on 3D terrains, comparing 
to flat agar plates yet sufficiently approximates man-made transport routes studied. 
We discuss how nutrients placed in destination sites affect performance of slime mould, 
how the mould navigates around elevations, and offer an automaton model of slime 

• mould which might explain variability of the protoplasmic routes constructed. We also 
^ ^ display results of scoping experiments with slime mould on 3D terrains of Russia and 
^ "H United Kingdom. 
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^ 1. Introduction 

VO 

^ The increase of long-distance travel and subsequent reconfiguration of vehicular and 

* social networks requires novel and unconventional approaches towards analysis of dy- 
^ namical processes in complex transport networks, routing and localisation of vehicular 

networks, optimisation of interactions between different parts of a transport network 
^ during scheduling of the road expansion, and shaping of transport network structure. 

Historical development and reconfiguration of major transport routes is somewhat anal- 
^ ogous to biological growth and pattern formation. Are man-made networks really similar 

^ to biological transport networks? Can living substrates develop better, in terms of min- 

imal distances and maximum coverage, transport networks than humans made? There 
is only one way to answer the questions: by undertaking laboratory experiments with 
a living substrate that produces transport networks. Ants and plant roots are first in- 
stances that come to mind, however it is somewhat difficult to experiment with these 
creatures. In 2001 Nakagaki with colleagues published a paper on shortest path approx- 
imation with live Plasmodium of Physarum polycephalum [17j. Three years later Tsuda, 
Aono and Gunji implemented Boolean logical gates with live slime mould [24j. Thus the 
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scientific world was introduced to a unique living computing substrate, which is easy to 
cultivate and handle, monitor and experiment with. 

Plasmodium is a vegetative stage of acellular slime mould P. polycephalum, a single cell 
with many nuclei, which feeds on microscopic particles [22j. When foraging for its food 
the Plasmodium propagates towards sources of food, surrounds them, secretes enzymes 
and digests the food; it may form a congregation of protoplasm covering the food source. 
When several sources of nutrients are scattered in the Plasmodium's range, the Plasmod- 
ium forms a network of protoplasmic tubes connecting the masses of protoplasm at the 
food sources. A structure of the protoplasmic networks is apparently optimal, in a sense 
that it covers all sources of nutrients and provides a robust and speedy transportation 
of nutrients and metabolites in the Plasmodium's body [TTt [T8t [T9]. 

To uncover analogies between biological and human-made transport networks and to 
project behavioural traits of biological networks onto development of vehicular trans- 
port networks we conducted a series of experimental laboratory studies on evaluation 
and approximation of motorway networks by P. polycephalum in fourteen geographi- 
cal regions: Africa, Australia, Belgium, Brazil, Canada, China, Germany, Iberia, Italy, 
Malaysia, Mexico, The Netherlands, UK, and USA [5j. We represented each region with 
an agar plate, imitated major urban areas with oat flakes, inoculated Plasmodium of P. 
polycephalum in a capital, and analysed structures of protoplasmic networks developed. 
For all regions studied [5], networks of protoplasmic tubes grown by Plasmodium match, 
at least partly, the networks of man-made transport routes. The shape of a country and 
the exact spatial distribution of urban areas, represented by sources of nutrients, may 
play a key role in determining the exact structure of the Plasmodium network. The ex- 
periments [5] had one obvious deficiency: slime mould propagated on a flat substrates, 
agar plates. Poor matching between protoplasmic networks and man-made transport 
networks in some countries could have been attributed to particular configurations of 
elevations and depressions in the countries's terrains. How does slime mould P. poly- 
cephalum behave on a 3D terrain? 

We envisaged that Plasmodium of P. polycephalum will not completely 'ignore' 3D 
topography because the slime mould is gravisensitive and positively geotropic ^ ^. 
Moreover, plasmodium shows morphological geopolarity, as discovered in [7]: the ec- 
toplasmic wall of a slime mould tube, lying or hanging on horizontal surface, is much 
thinner on side closer to earth. We expected that being placed on a 3D terrain with 
a source of nutrients slime mould would propagate towards the source of nutrients and 
navigate around elevations due to positive geotropism and relatively lower humidity of 
the elevations. 

Does slime mould on a 3D terrain match man-made transport routes better or worse 
than on a flat substrate? We conducted experiments on approximation of route 20 in 
USA and autobahn 7 in Germany with slime mould inoculated on a 3D plastic terrains. 
Why these specific routes? Route 20 (Fig. [l^) is the longest highway in USA, 5,415 km, 
it starts in Boston, Massachusetts, and end in Newport, Oregon (Fig.[l^) \2U\ [T4| [TBI [T2] . 
Route 20 crosses USA from east coast to west coast. Autobahn 7 (Fig. [l^) 963 km is 
the longest national motorway in Europe. It starts in Flensburg, passes via Hamburg, 
Hannover, Kassel, Fulda, Wiirzburg, Ulm, Kempten and end in Fiissen (Fig.[T]3) [9ll25]. 
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Figure 1. Schemes of route 20 (a) and autobahn 7 (b). (c) A protoplas- 
mic tube of P. polycephalum on a bare Nylon substrate. 



The paper is structured as follows. We described experimental setup in Sect. [2j 
Section [3] presents results of laboratory experiments of slime mould of 3D terrains and 
their analysis. Using excitable cellular automata models we imitate variability of slime 
mould routes in Sect. [Ij In Sect. [5] we test our approach to slime mould path based path 
finding on terrains of United Kingdom and Russia. A summary is provided in Sect. [6} 



2. Methods 



3D terrains of USA and Germany were ordered in http : //www . print ablegeography . 



[com/J where there were produced as follows. The elevation data are downloaded from 



DIVA-GIS (http://www.diva-gis.org/gdata), original source is CGIAR (http:// 



|srtm. csi.cgiar.org/). The data are projected with Mercator. Each terrain was 
printed using Selective Laser Sintered PA 2200 with Nylon 12. Terrain of Germany was 
8.2 cm wide, 11 cm long, and 3 cm high (dimensions measured between extremeties), 
and terrain of USA was 11 cm wide, 5.9 cm long and 3 cm high. 
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(c) 

Figure 2. Examples of slime mould routes developed on (a) USA shaped 
agar plate without attractants, (b) USA shaped agar plate with an oat 
flake placed in Boston, (c) 3D terrain of USA with an oat flake placed in 
Boston. 



In each experiment we inoculated Plasmodium of P. polycephalum in the positions 
of Newport, Oregon, USA (Fig. [l^) and Flensburg, Germany (Fig. [T]3), and oat flakes 
were placed in the positions of Boston, Massachusetts, USA, and Fiissen, Germany. 
The terrains were kept in closed yet naturally ventilated containers. The terrains were 
not artificially wetted or covered with any water retaining substrate but stayed in con- 
tainers filled with water (terrains were positioned above the water, without water even 
touching edges of the terrains). Humidity is proved to be sufficient for the slime mould 
to propagate on a bare nylon surface (Fig. [l]^). The terrains were kept in darkness, 
at temperature 22-27*^ C, except for observation and image recording. Configurations 
of Plasmodium networks were photographed with FujiFilm FinePix 6000 camera and 
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Figure 3. Examples of slime mould routes developed on (a) Germany 
shaped agar plate without attractants, (b) Germany shaped agar plate 
with an oat flake placed in Fiissen, (c) 3D terrain of Germany without 
attractants in Fiissen, and (d) 3D terrain of Germany with an oat flake 
placed in Fiissen. 
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scanned with an Epson Perfection 4490 scanner. Experimental setups with 3D terrains 
are illustrated in Figs. [Sj^d and[2]3. 

To understand influence of elevations on the slime mould's propagation we also con- 
ducted experiments on a flat agar gel. For these experiments we used 220 x 220 mm 
polystyrene square Petri dishes and 2% agar gel (Select agar, by Sigma Aldrich) as a 
substrate. Agar plates, about 2-3 mm in depth, were cut in shape of USA and Germany 
(Figs.|3^b and[2^b). 

We have conducted 11 experiments with USA-shaped flat agar, 14 experiments with 
3D terrain of USA, 13 experiments with Germany-shaped flat agar, and 16 experiments 
with 3D terrain of Germany. 

3. Results 

Usually it takes slime mould 4-5 days to cover the distance between Newport and 
Boston and 2-3 days to propagate from Flensburg to Fiissen (such difference in propa- 
gation time on the 3D terrains of almost the same size is possibly due to high amount of 
steep elevations in the West of USA, which prevent the slime mould from gaining speed 
and momentum when starting colonisation of the terrain). In experiments with Ger- 
many one day of slime mould's propagation roughly corresponds to 3-5 hours of real-life 
driving along the autobahn 7, and in experiments with USA to about 10 hours driving 
along route 20. 

Finding 1. Slime mould rarely reaches a destination which is not marked by a source 
of nutrients. 

In experiments without chemoattractants the Plasmodium was inoculated in start 
points, Flensburg and Newport, and no oat flakes were placed in destination sites, Fiissen 
and Boston. Due to proximal positions of start points the Plasmodium was only able 
to propagate south and east in USA (Fig. [2^), and south, east and west in Germany 
(Fig. [3^c). Schemes of the major protoplasmic tubes extracted from images of ten 
experiments on flat agar plates are shown in Fig. [3| The schemes demonstrate that 
slime mould's propagation was mainly limited by geometric shape of agar plate. There 
are clear evidences of the Plasmodium reflection from north and south boundaries of 
USA, and east and west boundaries of Germany. Such interaction with the boundaries 
of agar plates and randomised foraging behaviour of Plasmodium cause the slime mould 
to reach east cost of USA (Fig.|3^) and south boundary of Germany (Fig.[3]3) in arbitrary 
sites. Mishits — distances of Plasmodium's hits from destination sites divided by height 
of a country, in USA, or width in Germany — are tabulated in Fig. 

Finding 2. Slime mould approximates autobahn 7 in Germany 1.041 times better than 
it approximates route 20 in USA on agar plates. 

Slime mould's behaviour on flat agar plates with attractants, oat flakes placed in 
Boston and Fiissen was, rather partially, determined by chemo-attractants diffusing 
from oat flakes (Figs. [2]3 and[3]3). The Plasmodium demonstrated less foraging activity 
towards south in USA (Fig. [5^) and even less towards east and west in Germany (Fig.[5]3). 
Transport links built by slime mould between Newport and Boston are 1.046 times 
longer in average than a length of route 20 and the slime mould transport links between 
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(a) (b) 
Experiment Mishits USA Mishits Germany 



1 


0.187 


0.072 


2 


0.189 


0.094 


3 


0.276 


0.138 


4 


0.268 


0.156 


5 


0.426 


0.192 


6 


0.674 


0.249 


7 


0.726 


0.045 


8 


0.714 


0.091 


9 


0.659 


0.112 


10 


1.043 


0.509 


Mill 


0.187 


0.045 


Max 


1.043 


0.509 


Average 


0.516 


0.166 


Median 


0.543 


0.125 


Standard deviation 


0.289 

(c) 


0.135 



Figure 4. Results of experiments on growing shme mould without at- 
tractants, i.e. no oat flakes placed in destination sites Fiissen and Boston, 
(a) USA. (b) Germany. Major protoplasmic tubes extracted from dif- 
ferent experiments are shown by different colours (gradations of grey), 
(c) Statistics of mishits: deviations from destination sites divided by a 
heigh, USA, or width, Germany, of the country. 
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(c) (d) 

Figure 5. Results of experiments on approximation of route 20 in USA 
(ab) and autobahn 7 in Germany (cd) with shme mould. Slime mould is 
inoculated in Newport, Oregon (ab) and Flensburg (cd) and oat flakes are 
placed in Boston, Massachusetts (ab) and Fiissen (cd). Experiments are 
conducted on flat agar plates (ac) and 3D Nylon terrains (bd). Major 
protoplasmic tubes extracted from different experiments are shown by 
different colours (gradations of grey). 




Figure 6. Examples of slime mould navigating around mountains in 
USA (a) and Germany (b). 
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Table 1. Lengths of slime mould routes from start sites to destination 
sites divided by length of route 20 in USA and length of autobahn 7 in 
Germany. 



Experiment 


USA Agar 


USA 3D Germany Agar 


Germany 3D 


1 


0.963 


1.049 


1.005 


1.290 


2 


1.028 


1.056 


1.030 


1.117 


3 


0.924 


0.979 


0.962 


1.280 


4 


0.965 


1.278 


0.964 


1.010 


5 


0.954 


0.984 


1.016 


1.130 


6 


1.020 


0.956 


1.048 


1.252 


7 


1.009 


1.028 


1.016 


1.037 


8 


1.026 


1.108 


0.981 


1.065 


9 


1.147 


1.110 


0.994 


1.039 


10 


1.325 


1.156 


0.971 


1.093 


11 


1.148 


1.342 


1.014 


0.995 


12 




1.211 


1.093 


1.413 


13 




1.069 


0.967 


1.077 


14 




0.997 




1.040 


15 








1.359 


16 








1.327 


Min 


0.924 


0.956 


0.962 


0.995 


Max 


1.325 


1.342 


1.093 


1.413 


Average 


1.046 


1.095 


1.005 


1.158 


Median 


1.020 


1.063 


1.005 


1.105 


Standard deviation 


0.117 


0.116 


0.038 


0.139 



Flensburg and Fiissen are just 1.005 times longer in average than a length of autobahn 
7 (Tab. [T]). Assuming that chemo-attractants are not only airborne but also diffuse in 
a substrate, agar gel, shape of Germany provides favourable conditions for slime mould 
to detect diffusing chemoattractants: in USA the source of attractants is masked by 
indentations made by Lake Ontario, Quarry, Lake Huron and Lake Michigan (Fig. [5^). 
This may lead to prevalent diffusion of chemoattractants southwards and subsequent 
deviation of growing slime mould towards southern parts of USA. 
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Finding 3. Slime mould approximates route 20 in USA 1.095 times better than it ap- 
proximates autobahn 7 in Germany on 3D terrains. 

Table [T] shows statistics on lengths of major tubes representing route 20 and autobahn 
7 on 3D terrains. On average, slime mould routes are 1.095 times longer than route 20, 
and 1.158 longer than autobahn 7. In experiments with 3D terrains slime mould grows 
on bare Nylon terrains. Thus even if diffusion of chemoattractants occurred in the thin 
film of water covering the terrains, due to high humidity in closed containers, we believe 
the airborne attract ants would play a major role in guiding the Plasmodium toward 
destinations. Absence of any guidance force in a substrate together with non-flat terrain 
(Fig. [g]) cause the slime mould to deviate, sometimes substantially, from a shortest route 
towards its destination site (Figs. [2]3 and[3]l). 

In experiments with 3D terrain of USA slime mould is inoculated in Newport, Oregon. 
As soon as the Plasmodium recovers from inoculation procedure and begins exploration 
of a surrounding environment it finds that its way towards east is blocked by Cascade 
Range mountains. 

In six out of fourteen experiments with 3D terrain of USA slime mould (Fig. 
propagates north along Cascade Range mountains. If the 3D terrain included Canada 
the Plasmodium would propagate further north, however in present experiments the 
slime mould turns east along national boundary of USA. It crosses Rocky Mountains at 
the sites of Glacier National Park and Flathead National Forest. The Plasmodium then 
passes through North and South Dakota. In Minnesota plasmodia (recorded in different 
experiments) split: half of them continues into Wisconsin along interstate 94, another 
half turns south to Iowa, along interstate 45, and then turn east to Illinois along interstate 
70. On encountering Lake Michigan the 'Winsconsin bunch' of plasmodia also moves to 
Illinois. Scenario of further propagation is very common and registered in nine out of 
twelve experiments: the slime mould propagates through Indiana, Ohio, Pennsylvania, 
New York and, finally, enter the its destination site in Massachusetts (Fig. [s]^). 

In seven out of fourteen experiments with 3D terrain of USA slime mould propagates 
along bed of Columbia river and crosses Cascade Range between Gifford Pinchot and 
Mountain Hood National Forests (Fig. [5]3). The Plasmodium crosses Rocky Mountains 
range north of Yellowstone National Park, in the area between Butte and Bozeman. 
After traversing Rocky Mountains protoplasmic routes deviate: in different experiments 
slime mould can pass through Wyoming, Iowa, Minnesota. 

Only in one experiment we observed initially propagating along west coast of USA, 
west of Cascade Range, and along Sierra Nevada, parallel to interstate 5. The slime 
mould then turns east in region between Los Angeles and Las Vegas, move north-east, 
propagates along interstate 15. The slime mould passes nearby Salt Lake and in the 
region of Butte and Bozeman joins routes adopted by majority of slime moulds. 

On approaching Boston, most slime mould routes tightly match route 20, and only two 
protoplasmic tubes deviate east. One tube follows interstate 81 and another propagates 
along interstate 85. 

In contrast to USA, slime mould propagating on 3D terrain of Germany shows much 
higher variability in configurations of major protoplasmic routes (Fig. [5]l). On leaving 
Schleswig-Holstein slime moulds's trajectories dissipate west following autobahn 22 (five 
out of sixteen experiments) and east (ten our of sixteen experiments). Two of the 
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eastern routes propagate along autobahn 20, while others positions themselves between 
autobahn 7 and autobahn 26. Only in one experiment slime closely followed autobahn 
7 from the inoculation site. 

The western routes propagate along autobahns 36, 37 and 3, and one protoplasmic 
tube makes a detour south along autobahn 5 and then returning back north. Three 
out of five western routes, developed by slime mould, pass in the areas of Osnabriik, 
Miinster, Bonn, Mainz and join autobahn 7 in the region between Crailsheim and Hei- 
denheim. Two western routes join autobahn 7 in the region between Hildesheim and 
Kassel. Protoplasmic tubes spreading in eastern Germany join autobahn 7 either in the 
region of Wiizbourg or Fiissen (Fig. [5]l). 

Finding 4. Transport links developed by slime mould of P. polycephalum are longer 
than route 20 and autobahn 7. 

How often do slime mould build routes shorter than route 20 and autobahn 7. In 
Tab. [l] we find that slime mould propagating in USA builds transport links shorter than 
route 20 in circa 36% of experiments on flat agar and around 29% of experiments on 3D 
terrain. Slime mould growing on Germany-shaped flat agar develops protoplasmic tubes 
shorter than autobahn 7 in nearly 46% of experiments, while on 3D terrain in just one 
of fifteen, 7%, experiments. Possible explanations would be that road engineers were 
allowed to dig tunnels through mountains and built bridges while slime mould did not 
do that. 

4. Understanding slime mould via computer simulation 

Living systems are high variable, they rarely behave exactly the same way again and 
again even in rigid experimental setups. P. polycephalum does not make an exclusion: 
no two paths between start and destination sites are the same. Is it because of intrin- 
sic randomness of slime mould or for some other reasons? We believe exact routes of 
Plasmodium propagation are determined by an agility of Plasmodium, they depend on 
the Plasmodium's physiological state and level of metabolism. Active and vigour slime 
mould can climb up and cross elevations. Thus it undertakes less deviations and travels 
towards its destination site more close to as-crow-flies path. We show in computer sim- 
ulation that indeed topology of protoplasmic routes depends on a level of Plasmodium's 
activity. We represent 'activity' as a level of excitation in a two-dimensional cellular 
automaton. 

There are many ways to imitate growth and behaviour of P. polycephalum^ including 
flow networks [23j, numerical integration of partial differential equations [2J, multi-agent 
systems [llj and their combinations with cellular automata [lOl [2lJ. Based on our 
previous successful experience of imitating slime mould with localised wave-fronts in 
Oregonator model [2j we decided to imitate growth of the Plasmodium with excitation 
wave-fronts in cellular automaton and extract shortest path using pointers. Such ap- 
proach on approximation of a shortest path in cellular automata P works well not only 
in cellular automaton models but in crystallisation based computing devices [3J. 

We imitate active zones of growing slime mould in a 2D cellular automaton. The 
automaton is a rectangular array A of cells, finite state machines, which take four 
states: resting o, excited +, refractory — and precipitate #. Each cell x receives RGB 
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(c) e = 150 (d) e = 150 



Figure 7. (a) Time-lapse configuration of excitation wave-fronts prop- 
agating from Newport to Boston on a satellite image of USA. The wave- 
fronts are displayed at every 25tli step of iteration, (b) Configurations 
of pointers imposed by the wave-fronts in (a). Eight major orientations 
are represented by unique colours (gradations of grey). 6 — 150, d — 10, 
59 — 10. Cellular automaton has 500x253 cells. 

values {px^ l^x^lx) of the corresponding pixel x of a satellite image of the country under 
study (satellite images are taken in maps.google.com and reduced to sizes 500x253 
pixels, USA, and 505x707, Germany). Let be a state of cell x at time step t. All 
cells update their states simultaneously by the same rule, depending on states of their 
neighbours u{x) — {y ^ A \ \xy\ — 1}: 

i+, if = o and 5^ > and px < 0^ and 7^ < c and f3x < c 
, if = + 
#, ifx^ = - 
o, otherwise 

where 5^ is a sum of excited neighbours, 5^ = J2yeu{x)X{y^i~^)^ where xiu^i^) = 1 if 
= +, and otherwise; c = 230 and 9^ is updated by the rule 



9^ - 69, ifS^ <e 
9^ + 59, if 5^ > e + d 
9^, otherwise 
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where = X^^gA +)? ^ is a range of excitation activity. Value e represents a 
minimum level of activity or energy of imitated slime mould, and value e + d is a an 
upper limit of activity. A resting cell excites if it has at last one excited neighbour 
(5* > 0) and red value of a pixel (of satellite image) corresponding to cell x does not 
exceed 0^ and green and blue values do not exceed a constant c, the constant c is chosen 
230, so excitation wave does not propagate outside the country's shape (Fig. [T^c). The 
automaton behaves like a classical excitable automaton with the only difference that a 
refractory state is not changed back to a resting state but to an absorbing precipitate 
state. 

Value of 6^ (the same for all cells) is updated at very step t of iteration in such a 
manner to keep excitation propagating yet minimise propagation over high elevation 
sites (high elevations sites have higher values of px). Namely, when number of excited 
sites is less than a specified value e the threshold 0^ is increased to let excitation 
propagate over higher elevations and do not extinct. However, when there are too many 
excited sites, 5^ > e + d, the threshold 0^ is decreased and therefore excitation fronts 
can enter less elevations. 

Excitation wave-fronts imitate active growing zones of Plasmodium while pointers 
imitate protoplasmic pre-network. A thickest protoplasmic tube is selected by tracing a 
path along the pointers from destination site to start site. 

Initially all cells are resting but a cell corresponding to start site (Flensburg in Ger- 
many and Newport in USA) is excited. Excitation wave-front propagate from the start 
site to destination site avoiding elevations specified by 0^ (Fig.[7ic). Iterations are halted 
when either excitation reaches destination-site (Fiissen in Germany and Boston in USA) 
or becomes extinct. 

To trace a shortest path from a destination site to a source-site we must supply cells 
with unit vectors, or pointers to their neighbour, showing a local direction excitation 
came from. Each cell x is assigned a pointer p^; which has a zero length at the beginning 
of cellular automaton development. The pointer px becomes non-zero vector when cell 
X is excited: 



When we want to trace a shortest path — as presented by configuration of pointers 
developed in an an excitable medium — we start with destination site and then follow 
pointers to extract sites of the path one by one. In rare situations when path becomes 
loop locked we we choose next neighbour to move to in random and continue following 



Finding 5. A length of a path developed by the excitable cellular automaton decreases 
with increase of e. 

Value 0^ is increased when a level of excitation (number of excited cells) exceeds e + d. 
Thus the larger is e the less chances for 0^ to increase. Therefore, for large values of e 
more sites can be excited, and thus excitation wave fronts can climb on higher elevations. 
The dependency is illustrated in Fig. [9^. Lengths of paths calculated by the automaton 
drops substantially when e increases from 40 to 60, and then stabilises when e exceeds 




), ifx^ 



= + 



pointers (Figs, [s] and 10). 
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(i) e = 140 (j) e = 150 



Figure 8. Examples of paths between Newport area and Boston area 
developed by cellular automaton on a satellite image of USA. The follow- 
ing parameters are used: e is indicated in subfigures, d = 10, 0^ = 150 
and 66 = 10. Cellular automaton has 500x253 cells. 
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20 40 60 80 100 120 140 160 180 

(a) Route length 




Time 

(b) Number of excited sites 

Figure 9. Some integral characteristics of the model, illustrated in 
Fig. [sj (a) Lengths of routes approximated, (b) Dynamics of excita- 
tion level. The following parameters are used: d = 10, 9^ = 150, and 
69 = 10. 

100 (Fig.|9^). A path becomes shorter but at the costs of more energy spent on ascending 
elevations. 

Excitation propagates as a single wave- front only at the beginning of development, 
see high rise in a number of excited sites in Fig. [9]3. As soon as excitation level exceeds 
e + d value of 9^ decreases and thus less elevations can be climbed on. This leads to 
splitting of a single excitation wave-front to several fragments of wave-fronts whose size 
and configuration are controlled by a configuration of elevations in the fronts' vicinity. 
Thus for a low values of e we have many small wave-fronts scouting the satellite image 
(Fig. [9^). For high values of e we observe a rather classical wave, which sometimes gets 
split by elevations (Fig. ^ 
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Figure 10. Examples of paths between Flensburg area and Fiissen area 
developed by cellular automaton on a satellite image of Germany. The 
following parameters are used: 0^, e and d are indicated in the subfigures 
SO = 15. Cellular automaton has 505x707 cells. 



Figure [8] shows how topology of a path between Newport and Boston changes with 
increase of e from 30 to 150. For almost any path in (Fig. [s]) we can find a similar 
protoplasmic route obtained in experiments with slime mould (Fig. especially for 
routes lying along or south of route 20. Paths developed by propagating wave- fronts for 
high values of e reasonably close match route 20, see e.g. Fig. 

Dependence of a path length on activity level is demonstrated with cellular automaton 



developing path between Flensburg area and Fiissen area in Germany (Fig. 10). For 



low levels of activity, e.g. e = 20 (Fig. MOk) path propagates from Flensburg east to 
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Brandenburg, enters Sachsen, propagates west-south-west along Erzgebirge, reflected 
slightly north by Thiiringen Wald. The path navigate between Kellerwald, Vogelsberg, 
Taunus and Spessart, propagates south along Shwarzwald, and only on reaching south 
boundary of Germany turns towards its destination site. The path matches several 
plasmodial routes obtained in laboratory experiments (Fig. [5]l). 

When we increase e to 65 a path developed by excitable cellular automaton propagates 
from Schleswig-Holstein through Hannover to Bremen and Niedersachsen and navigates 
around Harz. The path turns east to Thiiringen, gets reflected by north part of Bohmer- 
Wald to south-west and then propagates to its destination-site (Fig. fiOjp). For high 
values of e path develops almost as a straight hne (Fig. [To}). 

In summary, plasmodium of P. polycephalum does not compute a shortest path per 
56, but the path optimal for the amount of resources involved. 



5. Scoping experiments witr 3D terrains of United Kingdom and Russia 
Example of one of the scoping experiments with 3D terrain of the United Kingdom 



is illustrated in Fig. 11 We used 3D terrain 6.2 cm wide, 11.0 cm long with maximum 
elevation 1.0 cm. We inoculated plasmodium in the northern part of Scotland, in the 
area between Tongue and Thurso, and placed an oat flake (which acts as a chemoat- 
tractant) to western part of England, in the area between Truro and Penzance (due to 
size of an average oat flake, c. 5 mm, it was impossible to place the flake in the location 
of the Land's End). A network of protoplasmic tubes developed does not match any 
major roads in the Highlands. Some matches between the roads and the tubes are ob- 
served when the slime mould approach Edinburgh in proximity of A84 and M9 (compare 



Fig. 11=1 and b). South of Edinburgh the slime mould propagates along Al and A697, 
passing close to Eyemouth. 

While propagating further south the plasmodium passes Durham, near Stockton- 
on-Tees, Leeds and Bradford area, Bolton and Crewe. On entering Wales the slime 
mould branches: one branch goes through Brecon to Carmarthen, another to the area of 
Cardiff and Newport (Fig.[TT|). The plasmodium crosses from Newport to Bristol, passes 
Bridgewater and then follows along M5. After passing Exeter the slime mould spreads 
all over Cornwall. Also, in the area between Coventry and Northampton a redundant 
branch appears. It follows towards location of London along Ml and partly A38 and 



M42, and then follows south between A3 and M3 (Fig. [IT]). 

Experimental image Fig. [TT| i shows well that while propagating from Scottish High- 
lands to Cornwall slime mould avoids entering or climbing on top of major elevations. 
Namely, the plasmodium navigates around Ben Nevis and Grampian Mountains, avoids 
Southern Uplands by passing close to Firth of Forth. It passes east of Cheviots Hills and 
crosses Pennines in the least elevated place, as if along a line connecting York to Manch- 
ester. The slime mould avoids Cotswold Hills by going from Newport area to Bristol. 
Its redundant branch passes west of Chiltern Hills while propagating to London. 

The experiment described above and illustrated in (Fig. [TT^ ) is quite an illustrative 
one. As you can see in the scheme (Fig. [TT] ^) of major protoplasmic tubes, extracted 
from images obtained in four experiments, the slime mould, in general, follows the same 
route as described above with minor variations in navigating around mountains. 
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Figure 11. Slime mould on 3D terrain of United Kingdom, (a) Scanned 
image of one experiment, (b) A scheme of major roads networks in 
UK (b) ^ 
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Figure 12. Images of four experimental setups, where slime mould de- 
velops a transport network between (ab) Moscow and Petropavlovsk- 
Kamchatsky and (cd) Moscow and Vladivostok, (e) Map of Trans- 
Siberian railway (red), Baikal- Amur mainline (green) [13j. 



ROUTE 20, AUTOBAHN 7, AND P. POLYCEPHALUM 21 




Figure 13. Example of slime mould navigating mountains, (a) West 
to east view of slime mould crossing Ural Mountains along a virtual line 
connecting Vorkuta to Labytnangi. (b) Slime mould propagating along 
slops of Kamchatka mountains and into Kamchatka peninsula, (c) South 
to north view of slime mould crossing Ural Mountains along a virtual 
line connecting Vorkuta to Labytnangi. (d) Slime mould passes through 
Central Siberian Plateau north of Enashimsky mountain, in the region of 
Tura city. 
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In experiments with Russia two terrains were used. One terrain, white, is 20.0 cm 
wide, 10.6 long and 5.0 cm high (at highest elevation, peaks of Elbrus were cut of), 
another, black in colour, has the same length and width yet, highest elevation is 1.5 
cm. Plasmodium was inoculated in a domain between Moscow and Kazan and oat flakes 
placed in locations of Petropavlovsk-Kamchatsky or Vladivostok (Fig. 12). 

Slime mould navigated towards sources of chemo-attractants often passing mountains 
in the places of lowest elevations (Fig. [131 ) demonstrated a high-degree of variability 
in developed protoplasmic transport links. For example, Plasmodium growing on a 
higher-elevated model (white terrain) tended to cross Ural Mountains in their northern 
parts (Fig. [T2]3d) while the Plasmodium growing on lower-elevated model (black terrain) 
avoided Ural Mountains in their southern parts (Fig. 12a.c). 

Only in one of the experiments illustrated the slime mould avoided Central Siberian 
Plateau by passing around its northern terminus (Fig. [T2|d) . In three other experiments 
the Plasmodium developed a protoplasmic link which either avoids the elevation by 
passing between Central Siberian Plateau and Sayan Mountains (Fig. \i2)p ) or crossing 
the Plateau along Podkamennaya Tunguska river (Fig. 12=ib). 

Few experiments undertaken with 3D terrain of Russia did not provide a substantial 
evidence that Plasmodium's transport networks connecting Moscow region and Vladi- 
vostok closely match Trans-Siberian railway or Baikal- Amur mainline (Fig. [12^). One 
experiment, shown in Fig. [T2|3, was a pleasant surprise: major protoplasmic tube grown 
along Trans-Siberian railway from Moscow-Kazan region to Omsk and then to Krasno- 
yarsk. The protoplasmic tube then propagates along Baikal- Amur mainline till turning 
south towards Vladivostok. 



6. Discussion 

Acellular slime mould Physarum polycephalum approximates shortest path in a labyrinth 
when pieces of Plasmodium are distributed in all channels of the labyrinth [17j . The slime 
mould also finds a shortest path in a labyrinth from a single start site while being assisted 
with diffusion gradient of chemo-attractants |4]. The slime mould satisfactory approx- 
imates man-made transport networks between major urban areas when these areas are 
represented by sources of nutrients [5j. All the above experiments were conducted on a 
flat agar gel. A few questions remained. How well the slime mould approximates a path 
on a substrate with elevations? Do slime mould paths grown on flat substrate differ from 
the paths developed on 3D models of countries? 

To find an answer we conducted a series of experimental laboratory studies on ap- 
proximation of route 20, the longest road in USA, and autobahn 7, the longest national 
motorway in Europe on 3D terrain plastic models of USA and Germany. 

We inoculated Plasmodium of P. polycephalum at one end-point of a road approxi- 
mated and placed sources of nutrients at another end-point. The slime mould grown 
from its inoculation site, explored 3D terrains, reached sources of nutrients and repre- 
sented a path between start site and destination site with its sickest protoplasmic tube. 
We have also done control experiments on a flat agar gel, and on substrates without 
sources of nutrients. 

In laboratory experiments we demonstrated that route 20 is approximated by live 
slime mould on 3D terrain closer than autobahn 7. Autobahn 7 is better approximated 
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on a flat agar gel. On flat surface average protoplasmic route in Germany is just 1.005 
times longer than autobahn 7 while an average protoplasmic tube in USA is just 1.046 
times longer than route 20. On 3D terrains of USA average slime route is 1.095 times 
longer than route 20, while on terrain of Germany average slime route is 1.158 times 
longer. It is not yet clear why plasmodium of P. polycephalum performs better on 
USA terrain. We can speculate this is due to location of transport routes approximates 
relative to national boundaries of the countries. Route 20 in USA leans towards north 
boundary of the country while autobahn 7 lies more or less equidistantly from east and 
west boundaries of Germany. Therefore in USA slime mould can substantially deviate 
from route 20 by propagating only south, in Germany growing slime mould can deviate 
west and east. 

In computer simulations we imitated active growing zones of plasmodium with excita- 
tion wave-fronts and structure of protoplasmic tubes with pointers in excitable cellular 
automata. We shown that conflguration of a path developed is determined by a level 
of activity of an imitated plasmodium, expressed via relation between minimum and 
maximum allowed levels of excitation and maximum heigh of elevations allowed to climb 
on. Higher levels of activity lead to shorter paths developed. 
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